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We develop a new framework of the deformed quasiparticle-random-phase approximation (QRPA) 
where the Skyrme density functional and the density-dependent pairing functional are consistently 
treated. Numerical applications are carried out for the isovector dipole and the isoscalar quadrupole 
modes in the spherical 20 O and in the deformed 26 Ne nuclei, and the effect of the momentum 
dependent terms of the Skyrme effective interaction for the energy-weighted sum rule is discussed. 
As a further application, we present for the first time the moments of inertia of 34 Mg and 36 Mg 
using the Thouless-Valatin procedure based on the self-consistent deformed QRPA, and show the 
applicability of our new calculation scheme not only for the vibrational modes but also for the 
rotational modes in neutron-rich nuclei. 

PACS numbers: 21.60.Jz; 21.10.Re; 21.60.Ev; 27.30. +t 
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I. INTRODUCTION 

Exploring nuclei far from the stability line is one of the 
most actively studied fields in nuclear physics. The ex- 
otic nuclei have revealed many features of atomic nuclei 
that are quite different from stable nuclei. Examples are 
the emergence of the neutron halo [lj and the skin Q 
structures, the soft dipole excitations [H, the modifica- 
tions of some magic numbers [l^EJ and the appearance 
of new magic numbers instead 6J, the onset of new re- 
gions of deformation Q ■ These new features are strongly 
connected with the presence of the loosely bound neu- 
trons and the coupling with the positive energy contin- 
uum states. Under the new environment, the nuclear 
many-body correlations such as the pairing and the de- 
formation could have also unique features in exotic nu- 
clei i, d ED, HI HI 0, HE Ell • 

For describing multipole responses in exotic nuclei 
in the region of medium-heavy systems, there have 
been many attempts employing the particle-hole ran- 
dom phase approximation (R PA) fl7l . 18. [l9 | or the 
^asiparticle-RPAfQRPA) 0, M, ^MMM, US HI 
Hi, Ha [H HI, HI Hi H|| on top of the self-consistent 
Hartree-Fock (HF) or Hartree-Fock-Bogoliubov (HFB) 
mean fields. (See Refs. [IE H3, HH for extensive lists 
of references concerning the self-consistent (Q)RPA and 
mean- field calculations.) These studies arc largely re- 
stricted to spherical systems. Recently, low-lying RPA 
modes in deformed neutron-rich nuclei have been inves- 
tigated by several groups [H EE EL El EI EI EE HI • 
These calculations, however, do not take into account the 
pairing correlations, or rely on the BCS approximation 
for pairing, which is inappropriate for describing the pair- 
ing correlations in drip line nuclei due to the unphysical 
nucleon gas problem [8[ . 

In order to discuss simultaneously effects of nuclear de- 



formation and pairing correlations including the unbound 
quasiparticle states, we have developed in Ref. [47j a cal- 
culation scheme that carried out the deformed QRPA cal- 
culation based on the coordinate-space HFB formalism. 
The residual interaction in the particle-hole (p-h) chan- 
nel was a simplified Skyrme interaction neglecting the 
momentum dependence 48[ while a deformed Woods- 
Saxon potential was employed for the mean field. In 
Ref. [43, one step further was accomplished by using 
a self-consistent Skyrme-HFB deformed mean field while 
the p-h residual interaction corresponding to the same 
Skyrme force was approximated by its Landau-Migdal 
(LM) form [H [5l|. However, the resulting energy- 
weighted sum rule (EWSR) for the isovector dipole re- 
sponse was not fulfilled very accurately. 

A full consistency between the QRPA and HFB cal- 
culations is required for a quantitative description of the 
multipole strengths in exotic nuclei. This means that 
the same effective interaction or the same energy den- 
sity functional is used for both calculations. We note 
that fully consistent HFB+QRPA calculations with the 
Gogny effective interaction for deformed nuclei are now 
available [HJ , but the use of a harmonic oscillator basis 
is a drawback for describing the unique spatial structure 
of quasiparticle wave functions near the Fermi level in 
neutron drip-line nuclei. 

We therefore develop in this article a new method 
for solving the Skyrmc-HFB-QRPA problem in deformed 
systems while keeping the full velocity dependence of the 
p-h residual interaction. The Skyrmc-HFB mean field 
is calculated in the coordinate-space representation. Nu- 
merical calculations are performed in order to investigate 
the effects of the explicit treatment of the momentum- 
dependent terms of the effective interaction on the EWSR 
of multipole responses and on the moments of inertia of 
deformed neutron-rich nuclei. The decoupling between 
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the spurious mode of translation and intrinsic excitations 
is reasonably obtained in these consistent calculations. 

The article is organized as follows: In Sec. UH the 
method is explained. In Sec. IIII1 we perform the numer- 
ical calculations and investigate properties of the isovec- 
tor/isoscalar dipole and the isoscalar quadrupole modes 
in 20 O, the isovector/isoscalar dipole modes in 26 Ne, and 
the moments of inertia of 34 Mg and 36 Mg as well as the 
isoscalar quadrupole mode. Sec. IIVI contains the conclu- 
sions. We summarize in the Appendix the main formulas 
relevant to the Skyrme energy density functional to show 
clearly what are the terms that arc included in our ap- 
proach. 



II. METHOD 

In the present section, we explain our method of the 
deformed QRPA based on the Skyrme density functional. 
We solve the HFB equations [H, |5^| 



the QRPA equation in the matrix formulation [5£ 



h q (r,cr) - X q 



h q {r,a) 
-(h q (r,a) - A«) 



(1) 



in coordinate space using cylindrical coordinates r = 
(p,z,4>). We assume axial and reflection symmetries. 
Here, q = v (neutron) or 7r (proton). For the mean- field 
Hamiltonian h, we employ the SkM* interaction [§4j in 
the present numerical applications. Details for express- 
ing the densities and currents in the cylindrical coordi- 
nate representation can be found in Refs. The 
pairing field is treated by using the density-dependent 
contact interaction 15711. 



Vpair i? •> T ) 



5{r-r'). (2) 



where Qo{r) denotes the isoscalar density and P a the 
spin exchange operator. Assuming time-reversal sym- 
metry and reflection symmetry with respect to the x — y 
plane, we have only to solve for positive O and positive 
z. We use the lattice mesh size Ap = Az = 0.6 fm and 
a box boundary condition at (/9 max = 9.9 fm, z max = 9.6 
fm) for 20 O and 26 Ne, and (/w = 9.9 fm, z max = 12 
fm) for Mg isotopes. The quasiparticle energy cutoff is 
chosen at E qp ^ cut = 60 MeV and the quasiparticle states 
up to f2 7r = 15/2 ± are included (for 26 Ne, we include 
the quasiparticle states up to fi 7r = 13/2 ± in order to 
compare with the results of Ref. \4^)- Our calculation 
scheme for solving the HFB equations is quite similar to 
Ref. [58j |. whereas the reflection symmetry was not im- 
posed in Ref. [5^ ]. 

Using the quasiparticle basis obtained as the self- 
consistent solution of the HFB equations |T]), we solve 



-4, 
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X a/3 
I aj3 
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■yS 

(3) 

The residual interaction in the particle-hole (p-h) channel 
appearing in the QRPA matrices A and B is derived from 
the Skyrme density functional, 



v ph (r,r') = 



"Sky 



Sg(r')Sg(r) ' 



(4) 



where we neglect the spin-orbit interaction term C t VJ in 
Eq. (|A4[) as well as the Coulomb interaction to reduce the 
computing time. Wc also drop the so-called "J 2 " term 
Cj in both the HFB and QRPA calculations. Thus, the 
residual interaction reads 

V ph (r, r') = (ao + a' r ■ r' + (b + b' r ■ r> • a')S(r - r') 
+ (oi + a\r ■ t' + (h + b\r ■ r')a ■ a') 
x (tf 2 S{r-r')+S{r-r')k 2 ) 
+ {a 2 + a' 2 r ■ t' + (b 2 + b' 2 r ■ r')a ■ a') 
X (fct • 5(r - r')k), (5) 

where fe = ("V - V')/2i and fet = -(V - V')/2i. The 
coefficients in Eq. ([5]) are given in Ref. [26[ (For simplicity, 
the coefficients ao,a' ,bo and b' here include the density 
dependent terms and rearrangement terms of — /g in 
Ref. pr|). We do not include the pairing rearrangement 
terms coming from the second derivative of the pairing 
functional £ pa ir with respect to the normal density g. 

On the other hand, the residual interaction in the 
particle-particle (p-p) channel is derived from the pairing 
functional constructed with the density-dependent con- 
tact interaction @, 



,(r,r') 



5g(r')5g(r) 



(6) 



This altogether coincides with Eq. ((2]). 

Because the full self-consistency between the static 
mean-field calculation and the dynamical calculation is 
broken by the above neglected terms, we renormalizc 
the residual interaction in the p-h channel by an over- 
all factor f p h to get the spurious K n = 0~ and 1~ modes 
(representing the center-of-mass motion), and = 1 + 
mode (representing the rotational motion in deformed 
nuclei) at zero energy (v p h — * fph ■ Vph)- We cut the 
two-quasiparticle (2qp) space at E a + Ep < 60 MeV due 
to the excessively demanding computer memory size and 
computing time for the model space consistent with that 
adopted in the HFB calculation. Accordingly, we need 
another factor f pp for the p-p channel. We determine 
this factor such that the spurious = + mode asso- 
ciated with the particle number fluctuation (represent- 
ing the pairing rotational mode) appears at zero energy 

[v pp ► fpp • Vpp ) . 
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FIG. 1: The 2qp cutoff-energy dependence of the B(E2 | 
) values (upper panel) and excitation energies E(2 + ) (lower 
panel) for the first excited 2 + state in 20 O. 
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III. RESULTS AND DISCUSSION 

A. 20 O 

The application of our new calculation scheme is pre- 
sented at first for a spherical system. In Ref. [1^ , detailed 
properties of 20 O were investigated using the continuum 
QRPA based on the Skyrme density functional. In the 
present subsection, we show our results for the isovector 
dipole and isoscalar quadrupole modes in 20 O following 
the discussions in Ref. [l^j. In the HFB calculations, the 
pairing strengths t' = -280 MeV-fm 3 and t' 3 = -18.75t' 
with 7 = 1 are employed as in Ref. [22j ■ With the choice 
7 = 1 in the pairing interaction, the pairing rearrange- 
ment terms vanish in the residual interaction. 

Because we use a larger mesh size and a smaller box, 
the obtained solution is not exactly the same as in 
Ref. [12]; the calculated total binding energy is 157.7 
MeV and the average neutron pairing gap is (A)„ = 1.92 
MeV. There are also differences in the QRPA calcula- 
tions between the two calculations: the boundary condi- 
tions, the 2qp cutoff energy and the treatment of the spin- 
dependent interactions (a ■ a' terms in Eq. ([5])). In the 
present calculation, the transition spin density is treated 
exactly. 

In Fig. Q] we show the 2qp cutoff energy dependence 
of the excitation energies and electric quadrupole transi- 
tion strengths of the first excited state. In this figure, we 
show the excitation energies of the = + , 1 + and 2 + 
states. If the spherical symmetry is preserved perfectly, 
these energies should be degenerate. In the actual calcu- 
lation, however, the spherical symmetry is broken due to 
the finite mesh size. Therefore, we can consider this dif- 
ference (~ 150 keV) as the numerical error caused by the 
discretization of the coordinates. The transition strength 
B(E2 1) is a sum of the intrinsic transition strengths to 
the K* = 0+,±l+ and ±2+ states. Around the cut- 



FIG. 2: Response functions for the isovector dipole (upper 
panel) and the isoscalar quadrupole (lower panel) operators 
in 20 O. The dotted and solid lines correspond to the Landau- 
Migdal approximation (LM) and the full calculation taken 
into account the momentum dependence explicitly (mom). 
The dashed line denotes the calculation without the a ■ a' in- 
teraction while keeping the momentum dependence explicitly. 
The transition strengths are smeared by a Lorentzian function 
with a width of T = 1 MeV. The arrow indicates the neutron 
emission threshold E t h = 9.19 MeV. 



off energy at 50 MeV, one can see that both E(2 + ) and 
B(E2 |) converge enough. 

We show in Fig.[2]the response functions for the isovec- 
tor (IV) dipole and the isoscalar (IS) quadrupole opera- 
tors 

N Z 7 N 

F % = e lE riyuf ft) - e-J2nY 1K (n), (7a) 



A 

A 



^!=E r i y «rft). ( 7b ) 

i 

and the corresponding response functions defined as 



= EE 



r/2 \(i\F. 



K 



7T (E - huj,) 2 + r 2 /4' 



(8) 



In this figure, we also show the response functions ob- 
tained by using the LM approximation. This approxi- 
mation treats only approximately the momentum depen- 
dence of the Skyrme p-h residual interaction. The re- 
sulting contact force is defined by the density-dependent 
Landau parameters Fq,Fq,Gq and Gq, which are de- 
termined by the parameters of the Skyrme effective in- 
teraction [5l|. The renormalization factors for the full 
QRPA calculation are f p h = 1.095 and f pp = 1.180, while 
fph = 0.815 in the LM approximation. 
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FIG. 3: Energy weighted sum of the IV dipole and IS 
quadrupole strength functions. The solid and dotted lines 
represent the calculation taking into account the momentum 
dependence explicitly (mom) and in the Landau-Migdal ap- 
proximation (LM), respectively. The horizontal lines show 
the classical Thomas- Reiche-Kuhn (TRK) , and the RPA sum 
rule values for the isovector case including the enhancement 
factor, mi = mf(l + «) (re = 0.32 in 20 O with the SkM* 
interaction ). For the isoscalar case, the EWSR corresponds 
to the classical sum-rule value. 
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FIG. 4: Response functions for the IS dipole operators 
Eqs. CEHJ) and (TTTJ in 20 O. 



hancement factor, mi = m\ (1 + k) 27j where k = 0.32, 
whereas the calculation in the LM approximation un- 
derestimates by 13% the EWSR value. For the IS 
quadrupole mode, the calculated sum satisfies 98.9% of 
the EWSR value, which is comparable to the values ob- 
tained in Ref. [22|]. The calculation in the LM approx- 
imation overestimates by 9.3% the EWSR value. This 
accuracy is as same as in Ref. [22l |. 

Next, we discuss the decoupling of the spurious state 
from the physical states. The IS dipole compression 
mode is sensitive to the admixture of the center-of-mass 
motion [6(| because the response function to the IS dipole 
operator 



For the IV dipole mode, the location of the giant res- 
onance is quite different between the calculation in the 
LM approximation and that taking fully into account the 
momentum dependence. The peak position is shifted up 
in the latter case. In Ref. [22|, the same tendency was 
obtained. However, the shape of the giant resonance is 
different between the two calculations. This difference 
comes from the a ■ a 1 terms of the p-h interaction which 
were omitted in Ref. Indeed, if we drop these terms 
in our calculation we obtain a two-peak structu re ( see 
Fig. ^ which is consistent with the result of Ref. [22| • 

For the IS quadrupole mode, we can see a prominent 
peak at 2 — 3 MeV, as well as the giant resonance at 
around 20 MeV. The low-lying state is sensitive to the 
momentum-dependent components of the force while the 
position of the giant resonance remains the same. In both 
calculations we obtain a small peak at 9 MeV just above 
the threshold, but the transition strength and the exci- 
tation energy are not affected. This result is consistent 
with Ref. [il. 

Figure [3] shows the partial sum of the energy weighted 
strength defined as 



W{E X ) 



E 

huji <E X 



M<^a t |o>P 



(9) 



For the IV dipole mode, the calculated sum up to 60 
MeV reaches 99.5% of the EWSR value including the en- 



Fl S K = J2r*Y 1K (r i ) 



(10) 



contains the strengths of both the spurious mode and 
physical intrinsic excitations. In order to see how much 
the spurious component mixes to the physical states, we 
compare the response function to the operator Eq. (fT0|) 
with that to the corrected operator [6ll.l62l] 



^IS(cor) 



(11) 



where 77 = ^(r ). If the spurious component is com- 
pletely removed from the intrinsic excitations, the calcu- 
lated response functions to the operators Eqs. (fTU|) and 
(fTTj) should coincide with each other. 

Figure [4] shows the response functions to the IS dipole 
operators Eqs. (fTTJ|) and (JTTJ). At around 17 MeV, we 
can see a slight difference between the two responses. In 
Ref. [26[ , the admixture of the spurious component was 
investigated in detail. Accurate results could be obtained 
by using a very large cutoff energy of 140 MeV in the lqp 
spectrum included in the fully self-consistent QRPA cal- 
culations. In the present work, the 2qp space is much 
smaller than in Ref. [26| and some of the residual inter- 
actions are not included. Though there is some room 
for further improvements, the results obtained here are 
rather satisfactory. 
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FIG. 5: Response functions for the isovector dipole operator 
in 26 Ne. The dotted, dashed and solid lines correspond to the 
K 71 = _ , K* — 1~ and total responses, respectively. The 
thin dotted line represents the response in the LM approxima- 
tion. For the K w = 1~ response, the transition strengths for 
the K 1 * = ±1~ states are summed up. The arrow indicates 
the neutron emission threshold E t h = 6.58 MeV. 
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In Ref. [4!| , we have studied the properties of the low- 
lying isovector resonance in deformed 26 Ne using the LM 
approximation. Although the strength function observed 
at RIKEN [H could be well reproduced, the EWSR was 
not satisfied very accurately. In the present subsection, 
we present the QRPA results where the momentum de- 
pendence of v p h is fully included and we show how the 
EWSR is fulfilled with the new method. The parameters 
in the HFB calculation are the same as in Ref. (49[, the 
difference being only in the treatment of the momentum 
dependent terms of interaction in the QRPA calculation. 

We show in Fig.[5]the response functions for the isovec- 
tor dipole mode. The renormalization factors for the 
QRPA calculation are f p h = 1.093 and f pp = 1.225, while 
fph = 0.919 in the LM approximation. Compared to the 
response functions obtained by using the LM approxi- 
mation, the excitation energies of both the low-lying and 
giant resonances are slightly shifted up while the overall 
structure remains similar. We can clearly see the res- 
onance structure below 10 MeV as experimentally ob- 
served [63]. The resonance is governed by the K n = 0~ 
state at 8.7 MeV and the = 1" states at 9.1 MeV and 
9.6 MeV. The microscopic structure of the A" 71 " = 0~ state 
is given in Fig. [5] and Table HI Here, the single-particle 
states are obtained by rcdiagonalizing the self-consistent 
single-particle Hamiltonian h of Eq. (fT|). As in the case 
of the LM approximation, this A' 71 " = 0~ state contains 
mainly the neutron lp-lh configuration (2si/2)~ 1 (2p 3 / 2 ), 
whose squared amplitude is 0.75. The K n = 1~ state 
at 9.1 MeV has also the same main component, with a 
weight of 0.83. The state at 9.6 MeV is mainly generated 
by the (2s 1 / 2 ) _1 (2pi/2) excitation with a weight of 0.90. 
A difference between the results of the LM approximation 



FIG. 6: Two-quasiparticle (particle-hole) excitations gener- 
ating the K* = 0" state at 8.68 MeV in 26 Ne. The single- 
particle levels are labeled with the asymptotic quantum num- 
bers [NnsA]fl. The solid and dotted lines stand for the posi- 
tive and negative parities, respectively. The chemical poten- 
tial A is indicated by the two-dotted line. 
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FIG. 7: Same as Fig. [3] but for 26 Ne. 



and the present results is that the transition strength to 
the K n = 1~ state at 9.6 MeV (0.08e 2 fm 2 ) is larger than 
that to the state at 9.1 MeV (0.04e 2 fm 2 ). 

Figure [7] shows the energy weighted sum of the isovec- 
tor dipole strength function together with the sum rule 
values represented by the horizontal lines. The energy- 
weighted sum up to 60 MeV overestimates by only 1.6% 
the EWSR value including the enhancement factor n = 
0.32. In the calculation with the LM approximation, 
the energy-weighted sum is underestimated by about 
10% [4!|. This suggests that treating the momentum 
dependence of the Skyrme force explicitly in the QRPA 
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TABLE I: QRPA amplitudes for the = CT state at 



deformed systems 



8.68 MeV in 
B(Q IV 1) 



3 Ne. 



This mode has the isovector strength 
sum of backward-going 
Only components with 



8.72 x 10" 2 e 2 fm 2 , and the sum of backward-going 



amplitudes ]T \Y a/3 \ 2 = 5.24 x 10 

X^p — Y^p > 0.001 are listed. In the rows (h) and (i), the label 
z/1/2 - denotes a non-resonant discretized continuum state of 
neutron = l/2~ level. 
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FIG. 8: Same as Fig. g] but for 26 Ne. 
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calculation is crucial for satisfying the EWSR in the de- 
formed systems as in the spherical systems [22| ■ Because 
of the severe cutoff in the 2qp excitation energy, we can- 
not describe properly the energy region higher than the 
giant resonance. It would be interesting to see if calcu- 
lations in a larger 2qp space including the residual spin- 
orbit and the Coulomb interactions can improve the over- 
shooting of the EWSR value. 

Next, we discuss the IS dipole response. The corrected 
IS dipole operator Eq. (fTTj) with r\ = f^ 2 } is valid only 
for the spherical systems. We can extend it for the de- 
formed systems by following the discussion in Appendix 
of Rcf . [62| ■ One obtains the correction factor r\ for the 




(12) 



These correction factors coincide with rj = ^(r 2 ) in the 
spherical limit. 

In Fig. we show the response functions to the IS 
dipole operators with/without the correction. For the 
lowest state at 6.5 MeV, we can see a difference between 
the two calculations. However, the overall structures are 
not very different. We can consider that the spurious 
component is well removed from the pygmy resonance 
and the giant resonance. 



'Mg 



We show in Fig. [9] the response functions for the IS 
quadrupole transition in 24 Mg. We employ the same 
effective interactions for the HFB+QRPA calculation 
as in the 20 O nucleus. The rcnormalization factor is 
fph = 1-164. The giant resonance appears at around 
15-25 MeV. Since the ground state is prolately de- 
formed in our calculation ((3-2 = 0.4), we can see a clear 
if-splitting. Below 5 MeV, we can see a prominent peak 
for the = 2 + excitation. These are consistent with 
the fully self-consistent deformed QRPA calculation us- 
ing the Gogny force [52j. 

In Fig. 1101 we show the low-lying excitation spectrum. 
Here excitation energies are evaluated by [68[ 



E(I,K) = ^ RPA 



2Jtv 



(I(I+l)-K% (13) 



in terms of the vibrational frequencies, wrpa, and the 
Thouless-Valatin moment of inertia [65[ , Jrv , calculated 
microscopically by the QRPA. The calculated moment of 
inertia is large compared to the experimental rotational 
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FIG. 9: Response functions for the IS quadrupole transition 
in 24 Mg for the K* = + , 1 + and 2 + excitations represented 
by the solid, dotted and dashed lines, respectively. 
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FIG. 10: Excitation energy spectrum obtained by the QRPA 
calculation and the available experimental data [691 ]. 



band. This is because the present pairing interaction 
leads to vanishing neutron and proton pairing gaps. Fur- 
thermore, due to the absence of the coupling mechanism 
between the (3 vibration and the pairing vibration, the 
excited = + mode cannot acquire a substantial col- 
lectivity and the excitation energy remains large. The 7 
vibrational mode reasonably reproduces the experiment. 
This K T = 2 + state is mainly generated by the neutron 
and proton p-h excitations [211]3/2 -> [211] 1/2. Their 
contributions are 42% and 56%, respectively and the to- 
tal transition strength is exhausted about 55% by the 
two p-h excitations. The rest of the transition strength 
comes from the coupling to the giant resonance. 

D. 34 Mg and 36 Mg 

In Ref. (64|, the properties of the low-lying = + 
mode in 34 Mg have been studied and the moments of in- 
ertia were calculated using the Thouless-Valatin method. 
In the present subsection, effects of the momentum de- 
pendent components of the Skyrme interaction on the 



TABLE II: Ground state properties of 34 Mg and 36 Mg ob- 
tained by the deformed HFB calculation with the SkM* in- 
teraction and the mixed- type pairing interaction. Chemical 
potentials, deformations, average pairing gaps, root-mean- 
square radii for neutrons and protons, and the total binding 
energies are listed. 
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FIG. 11: Response functions for the IS quadrupole operator 
in 34 Mg (upper panel) and 36 Mg (lower panel). The arrows 
indicate the neutron emission thresholds E t h = 6.13 MeV and 
5.18 MeV for 34 Mg and 36 Mg, respectively. 



isoscalar quadrupole strengths and the moments of iner- 
tia are discussed. 

In the HFB calculations, the pairing strength param- 
eter is determined so as to reproduce the experimental 
pairing gap of 34 Mg (Agxn = 1.7 MeV) obtained by 
the three-point formula [66|. In Table [Til the ground 
state properties of 34 > 36 Mg are summarized. The strength 
4q = — 295MeV-fm 3 for the mixed-type interaction with 
7=1 leads to the pairing gaps (A„) = 1.71 MeV in 34 Mg 
and 36 Mg. We obtain for the proton intrinsic quadrupole 
moments Qo the values 62.2e-fm 2 and 60.1e-fm 2 in 34 Mg 
and 36 Mg, respectively. The reduced transition proba- 
bilities [63 are then B(E2;0+ -> 2+) = 5/16tt • Q 2 , = 
385e 2 -fm 4 and 359 e 2 -fm 4 in 34 Mg and 36 Mg. In 34 Mg, 
the neutron occupation probability of the [202]3/2 level 
coming up from the 1^3/2 orbit is 0.28, while that of the 
[321]3/2 level coming down from the I/7/2 orbit is 0.67. 
This approximately corresponds to the (l/7/2) 4 (1^5/2)~ 2 
configuration in the spherical shell model language. In 
36 Mg, the occupation probability of the ^[202]3/2 level 
becomes 0.64. These probabilities are consistent with the 
shell model results [70j . 

In the QRPA calculations of 34 Mg and 36 Mg, we cut 
the 2qp space at 60 MeV as in the previous subsections. 
We have checked that the transition strength to the ex- 
cited + state and its excitation energy converge at 50 
MeV cutoff. In the present case, the dimension of the 
QRPA matrix © for the = 0+ channel in 36 Mg is 
about 15,000, the memory size is 20.8 GB, and the CPU 
time is about 154,000s per iteration for determining the 
renormalization factor f pp using the SX-8 supercomputer 
at the RCNP. 

Figure 1111 shows the response functions for the IS 
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FIG. 12: QRPA strength distributions for (a) the isoscalar 
K n — + quadrupole p-h excitations and (b) the monopole- 
and (c) the quadrupole-pair excitations in 34 Mg. For compar- 
ison, unperturbed 2qp transition strengths are shown in the 
lower panels. 



quadrupole operator in 34 Mg and 36 Mg. The renormal- 
ization factors for 34 Mg and 36 Mg are f p h = 1.146 and 
1.139, respectively. In the pp channel, the renormaliza- 
tion factors are f pp = 1.197 and 1.215 for 34 Mg and 36 Mg. 
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FIG. 13: Excitation energy spectrum for the positive par- 
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ity states in 34 Mg and 36 Mg, and the available experimental 
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P20 = I drr 2 Y 20 (f)i>Ur T )#(r |). 
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In both nuclei, we can see a low-lying peak at 2 MeV 
and a three-peak giant resonance at 15 — 20 MeV. This 
three-peak structure corresponds to the giant resonance 
for the = + , 1 + and 2 + excitations. Because the 
deformation of 34 Mg is larger than that of 36 Mg, the K 
splitting is larger in 34 Mg. On the same figure are shown 
for comparison the results of the LM approximation. As 
in 20 O, the low- lying state is sensitive to the treatment of 
the momentum-dependent interactions while the position 
of the giant resonance is not much affected. 

The calculated energy-weighted sums up to 60 MeV 
for the K* = 0+ excitation in 34 Mg and 36 Mg amount 
to 99.0%, whereas in the LM approximation, they are 
overestimated by 11.5% and 11.0% the EWSR in 34 Mg 
and 36 Mg. This confirms that the EWSR is well satisfied 
in the present calculation. 

In Ref. [64| , we have also discussed the generic feature 
of the low-lying if* = + modes in deformed neutron- 
rich nuclei: In a deformed system where the up-sloping 
oblate-type and the down-sloping prolate-type orbitals 
exist near the Fermi level, one obtains a low-lying mode 
possessing enhanced strengths both for the quadrupole p- 
h transition and for the quadrupole p-p (pair) transition 
induced by the pairing fluctuations. In Fig. 1121 we show 
the strength distributions in 34 Mg of the quadrupole p- 
h, the monopole p-p and the quadrupole p-p transitions 



At 2.65 MeV, we obtain the collective K v = + mode 
possessing about 30 Weisskopf units for the intrinsic 
isoscalar quadrupole transition strength. (The electric 
transition strength is B(E2; 0^ -> 2f) = 13.0e 2 fm 4 .) 
The transition strength is enhanced by 10.6 times as com- 
pared to the unperturbed transition strength. For the 
quadrupole pair transition, the strength to this collec- 
tive state is enhanced by 14.9 times with respect to the 
unperturbed one, while the strength is not changed for 
the monopole pair transition. We have checked that the 
low- lying if" = + mode at 2.65 MeV is well decoupled 
from the pairing rotation. The transition strength for the 
number operator to this state is |(A|iV|0)| 2 = 2.87 x 10~ 5 . 

In Fig. [13l we show the low-lying excitation spectrum 
below 4 MeV for the positive parity states in 34 Mg and 
36 Mg. The Thouless-Valatin moment of inertia of 34 Mg 
is Jtv/Ti 2 = 4.23 MeV -1 when the LM approxima- 
tion is used, and 4.26 MeV -1 when the full momentum- 
dependent interaction is treated in the QRPA, while 
the Inglis-Bclyacv moment of inertia is jTeciyacv — 3.89 
MeV -1 . Due to the time odd components in the residual 
interactions ([5]) and ([6]), the moment of inertia Jtv be- 
comes about 10 % larger than Jneiyaev In 36 Mg, we ob- 
tain Jrv/h 2 = 4.20 MeV~\ and 4.24 MeV" 1 in the LM 
approximation. If we turn off the residual interactions 
we obtain 3.84 MeV" 1 . For both nuclei, the moments 
of inertia calculated by using the LM approximation are 
close to the results of the QRPA with the full velocity- 
dependent force. 
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IV. CONCLUSIONS 

We have developed a new calculation scheme of the de- 
formed QRPA using the Skyrmc density functional. This 
scheme allows one to include in the p-h residual inter- 
action the full velocity dependence of the Skyrme force. 
The only components of v p h not treated here are the spin- 
orbit and Coulomb two-body forces. Numerical applica- 
tions have been performed for some spherical and de- 
formed neutron-rich nuclei employing the Skyrme SkM* 
density functional and a local pairing functional for gen- 
erating the HFB mean field, pairing field and the resid- 
ual interactions in the QRPA calculations. In both the 
spherical and the deformed cases, we have checked that 
the energy-weighted sum rules for the isoscalar and the 
isovector operators are well satisfied. There is a distinct 
improvement over the sum rules predicted by the LM 
approximation, even in the case of isoscalar excitations, 
thus indicating the importance of full self-consistency in 
HF(B)-(Q)RPA calculations. 

Thus, this method enables one to describe multipole 
strengths quantitatively in deformed nuclei located in a 
wide range of the nuclear chart, even near drip lines. It 
has been also shown that one can apply it not only to 
the vibrational modes but also to rotational modes by 
employing the Thouless-Valatin procedure. 

Methods for solving the HFB-QRPA problem in de- 
formed systems using the Skyrme plus local pairing den- 
sity functional in a fully consistent way are still scarce. 
We have proposed here a new method and we have 
demonstrated its feasibility on some examples. This 
method has some advantages and drawbacks. One main 
advantage is the choice of solving the deformed HFB 
problem on a grid in coordinate space. This avoids ex- 
panding quasiparticle wave functions on a harmonic oscil- 
lator basis and introducing inaccuracies inherent to ex- 
pansions of loosely bound, or unbound wave functions 
on such basis. This may be of some importance when 
studying near drip-line nuclei. A practical drawback is 
the necessity of using a relatively large 2qp cutoff, and 
therefore computing times and memory storage are high. 
Our numerical studies show that a rather good accuracy 
is already reached if the 2qp energy cutoff is set at 60 
MeV. Doubtless that, in the future the capacity of com- 
puting facilities will largely improve and the 2qp space 
can be easily enlarged. 
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APPENDIX A: THE SKYRME DENSITY 
FUNCTIONAL 

The total energy of the system consists of the kinetic 
energy £kin, the Skyrme interaction energy £sky, the 
Coulomb energy £coui, the pairing energy £ pa i r and the 
correction of center-of-mass motion and rotational mo- 
tion £ corr ; 

£ = £kin + ^Sky + ^Coul + £pair + £corr, (Al) 

The kinetic energy is given by 

fldn = J ^^(r), (A2) 

where r is the kinetic density. In the present paper, we 
perform the numerical calculations using the SkM* in- 
teraction [54[, so the center-of-mass correction is just to 
replace the nucleon mass 1/m — > 1/m x (1 — 1/A). The 
correction for the rotational motion is not taken into ac- 
count. 

The Skyrme interaction energy is given as (z!, [74[ 

£sk y = / drH S ky(r), (A3) 



n Sky (r) = {C? Mr)k 2 (r) + Cf [ft,(r)]*?(r) 
t=o,i 

+ C t Ap g t (r)Ag t (r) + C As s t (r) ■ As t (r) 
+ CUgt(r)T t (r)-jf(r)) 
+ C?(s t (r).T t (r)-y 2 t (r)) 

+ C? J (g t (r)V ■ J t (r) + s t (r) ■ V x j t (r))}, 

(A4) 

where g denotes the nucleon density, s the spin density, T 

< — > 

the kinetic spin density, j the current tensor, J the spin- 
current tensor and J the spin-orbit current. All densities 
arc labelled by an isospin index t where t is (isoscalar) 
or 1 (isovector) and we assume no isospin mixing. 
The Coulomb energy is given as 

£coui = J drHcou\(r), 

H Coul (r) = dr'Mr)^ - (i) V( r)i 

(A5) 

where the exchange term in the Coulomb energy is 
treated in the Slater approximation 17511, and the higher 
order correction is found to be small [76|]. We follow the 
procedure of Ref. [77} for calculating the Coulomb poten- 
tial. 
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When wc use for the pairing interaction the form of 
Eq. ([2]), the pairing energy is given as 



where g denotes the abnormal (pairing) density and s 
the spin pairing density. 



f . — 



drH pair (r), 



H 



■pair 



■(r) = I 



£(e?(r)-S?(r)), (A6) 
t=o,i 
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